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This study presents an analysis of the fixed biases and variable 
errors affecting gun fire. The emphasis is placed on automatic cannons 
being developed for use as main guns for armored vehicles. 

It is proposed that barrel bend has been overlooked as a primary 
cause of inaccurate fire and that automatic cannons will be particularly 
susceptible to this condition. A formula is derived for approximating 
the magnitude of bend that can occur. 

A model is presented for estimating miss distance in the plane of 
the target as a result of fixed biases (including barrel bend) and other 
variable errors. A parametric analysis demonstrates the magnitude of 
resulting miss distance. Methods of avoiding the effects of barrel 
bend are presented. 
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I . INTRODUCTION 



With the introduction of vastly improved threat weapons and observed 
numerical imbalances, the Department of Defense over the last decade 
initiated many programs to analyze existing and advanced state of the 
art technologies for armored combat vehicles of the future . According 
to the United States Army Armor and Engineer Board E3. these programs 
were designed to improve the firepower, mobility/agility and protection 
of future combat vehicles while concurrently reducing weight and cost. 

As noted, an integral feature of these technological programs is the 
development of technology for improved firepower and current interest is 
vested in the medium caliber automatic gun. The Department of Defense 
analysts, especially research and development analysts, as well as 
tactical and strategic planners must evaluate data on the weapon per- 
formance. This need for methods of evaluating weapons systems becomes 
paramount since our future depends on intelligent and informed decisions 
at the earliest stages of any program. 

A current example began early in the 1970 ' s as the Defense Advance 
Research Projects Agency (DARPA) developed a program to demonstrate new 
technologies for future tanks. One of the major component technology 
areas in the program is a vehicle characterized by high mobility and 
agility and hence the program named, HIMAG. A critical element in the 
HIMAG program is the demonstration of the feasibility of a medium cali- 
ber automatic gun. The weapon pod proposed for evaluation is one 
containing the ARES 75 mm, solid propellant, rapid fire cannon using 
APFSDS and HE inert rounds. This 75mm gun is designed to fire in 
various modes and incorporates the capability for both rapid or single 
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shot fire. Here we have an explicit need for a method to evaluate the 
proposed weapon system in a burst mode to demonstrate the technical 
feasibility of the advanced automatic cannon. 

However, Pfleger and Cerrato H report that to date there has been 
only limited effort directed towards providing reliable, accurate anal- 
ysis of the performance of automatic gun systems. Primary emphasis has 
traditionally been to evaluate single shot weapon systems. Consequently, 
with the advent of the automatic cannon and the current need for effec- 
tive weapon analysis as generated by the Department of Defense , this 
paper's objectives are to: 

1. Develop a mathematical model to predict approximate miss distance _ 
of the projectile at the target plane for medium caliber automatic 
cannons in order to provide accurate and reliable information upon 
which hit probabilities can be developed for analysis by DOD. 

2. Develop a mathematical model to approximate barrel bend and use this 
model to determine the possible amount of error caused by barrel bend 
using the ARPA 75mm automatic cannon as a test case . 

3. Use the mathematical models developed in objectives "1" and "2" in 
a computerized parametric analysis of the AEPA 75nim cannon under a 
baseline M60A1 fire control system approximation and the more advanced 
XM-1 fire control system approximation. 

A more detailed statement and presentation of the automatic cannon 
accuracy evaluation problems are discussed in Chapter II, Nature of the 
Problem. The body of the study is presented in Chapter III, Presenta- 
tion and Analysis of Data. Chapter III, Section A, Miss Distance 
Analysis, provides a background for the reader on error sources in- 
volved with automatic cannons on combat vehicles as well as develops 
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the mathematical model on miss distance. Chapter III, Section B, 

Barrel Bend Analysis, provides a detailed analysis and mathematical 
model on barrel bend approximation. Once these mathematical models 
have been formulated, then Chapter IV, The Simulation Model, employs 
the models developed using the ARPA 75 mm automatic cannon as a test 
case under several assumed weather conditions. Lastly, Chapter V is 
dedicated to a summary and conclusions of the study as generated during 
the research effort and from analysis of the simulation. 
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II. NATURE OF THE PROBLEM 



Nearly 3000 tanks were destroyed during the 1973 Arab-Israeli war. 
Post-war analysis has shown that the clear majority of these tanks were 
destroyed by the main gun of other tanks. Thus, as late as 1973 » the 
tank remains the primary anti-tank weapon. Analysis also reveals the 
total rounds fired considerably exceeds the quantity expected to be 
fired for the number of hits actually achieved. It can be inferred 
from the data that tank guns are not nearly so accurate as commonly 
believed . 

In the event of a conventional war, armored vehicles constitute a 
primary threat to the security of our allies. If the tank is to be the 
primary anti-tank weapon in the future or if the tank is even to re- 
main effective as an anti-tank weapon, something must be done to im- 
prove the accuracy of tank gunfire . 

The primary effort toward improved accuracy has been the development 
of increasingly complex and sophisticated fire control instrumentation. 
Significant progress has been made in minimizing the effect of some of 
the factors which contribute to miss distance at the target, e.g. range 
estimation . 

A second approach to improving the probability of hit is currently 
being investigated. That is, to fire a burst of several rounds at the 
target instead of a single round. If each round of a burst can be fired 
in precisely the same intended trajectory, round-to-round dispersion 
will cause a shotgun effect at the plane of the target. The probability 
of a hit in such a burst would clearly be higher than for a single 
round . 
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To demonstrate this approach a weapon must be developed that is 
capable of burst fire. This means firing anti-tank rounds at rates of 
fire approaching 120 rounds per minute or, more appropriately, 2 rounds 
per second. Weapons of this sort are being developed. The ARES 75mm 
gun currently in development has similar characteristics. This gun will 
be used as a model when, in the course of this paper, it is necessary 
to utilize the specific characteristics of a rapid fire gun. 

Given a gun with rapid or burst fire capability, it must still be 
demonstrated that each round of a burst is fired with the same degree of 
accuracy. This requires the barrel to be pointing in the same direction. 
It must be demonstrated that the recoil or trunnion forces can be trans- 
ferred through the entire vehicle and the gun platform and gun can be 
returned to the proper position in the time available. Rapid fire also 
means higher barrel temperatures. High temperatures can affect barrel 
life and accuracy. Safety is another consideration in view of the 
possibility of ignition of the propellant due to chamber temperature or 
"cook offs." It is clear that there is much to be considered and 
demonstrated in the development of a rapid fire tank gun. 

A third approach is implied in the foregoing but should be addressed 
separately. Determine the factors which contribute to miss distance 
that have not been eliminated or minimized to date and find ways to 
compensate for or avoid their effects. It is felt that far too little 
attention is being given to two factors which make significant con- 
tribution to miss distance, i.e. round -to- round dispersion and barrel 
bend . 

Round- to-round dispersion is primarily the result of differences in 
the physical characteristics of each round which affect the trajectory 
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of the round. This dispersion is well known to those familiar with 
gunfire in general. Unfortunately, this error has been deemed "random" 
which implies "unavoidable." It is felt that further efforts to mini- 
mize this factor are justified. 

Barrel bend is a little known factor. Few recognize that gun 
barrels do, in fact, bend. Even fewer have considered the possible 
magnitude of the resulting miss distance. It is understood that weapons 
are being developed by West Germany and Great Britain that may avoid 
barrel bend or its effect. There is no information available to 
indicate that such a weapon is being considered in this country. 

A primary purpose of this paper is to demonstrate that barrel bend 
can cause a significant error, particularly in view of the higher 
barrel temperatures that will occur with rapid fire systems. 
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III. PRESENTATION AND ANALYSIS OF DATA 



A. MISS DISTANCE ANALYSIS 
1 . Introduction 

This section describes an analytical method formulated for use 
in the evaluation of the performance of automatic gun systems for combat 
vehicles against stationary ground targets. The need for such analytical 
effort is readily apparent, especially in view of the current emphasis 
on armored combat vehicle technology and the HIMAG program evaluating 
the use of the 75 ^ automatic gun . 

The first step in the analysis was to determine a significant 
effectiveness measure for the automatic gun. Traditionally, hit proba- 
bility (P^) has been used in this case and considerable effort in 
n 

measuring this effectiveness has been accomplished and guidelines 
established for first round or single shot hit probability by Frankford 
Arsenal £2,4]. However, with the advent of the automatic gun which can 
fire in a burst or rapid fire mode, the hit probability of subsequent 
rounds fired in rapid succession must be addressed. 

Secondly, it is felt that any analysis in this area must begin 
by determining the impact point of each round fired, either by single 
shot or burst fire. By mathematically modeling the impact point of each 
round fired at the target plane in relation to the center of the target,- 
a base would be provided for any additional analysis in hit probability 
evaluations. Logically, the farther the impact point of a round from 
the center of target, the lower the overall hit probability. The use 
of all error sources is designed to simulate the conditions which can 
be expected under a variety of actual combat conditions . According to 
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Brodkin DO , this measure, Quasi-Combat Condition, is intended to reflect 
as realistically as possible the combat potential of the system in re- 
gard to delivery accuracy. 

Lastly, under the scenario input to the simulation in Chapter IV, 
two terms mentioned by Cerrato 0 necessitate reiteration. An "occasion" 
is one tactical engagement between vehicles or stationary implacement. 

A "burst" is a series of rounds fired in rapid succession with the same 
point of aim. 

2 . Classification of Error Types 

Since the standard target is two-dimensional, the effects of 
both horizontal (X) and vertical (Y) error sources must be considered. 
These gunnery error sources cause certain errors in different magnitudes 
and in both directions at the target plane . Pfleger and 3ibbero DO re- 
port that this leads to the assumption that the horizontal and vertical 
errors are independent of each other. This method is valid if the 
target is aligned with the coordinates of the plane in which the errors 
are considered. 

From the classification of error types presented in Frankford 
Arsenal [^>5^ » gunnery errors can be classified into three types with 
respect to automatic cannons on combat vehicles . 

a. Fixed Biases are a consequence of system design and use, and 
are constant for a given range. In the fire control systems in use 
today, correction is made for many of these fixed biases. 

b. Variable Biases are considered to be constant for a specific 
engagement and/or burst of fire. They vary in a random fashion but not 
from round to round. Included in this category are: 

(i) Occasion-to-Occasion biases are errors which change over 
time, but change is so slow that they can be considered a constant value 
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over the period of an engagement . Generally, these errors are induced 
by the weapon system design and the tactical situation. 

(2) Burst- to-Burst biases are considered to have different 
values for each burst fired dinring the engagement and are considered to 
vary in a random manner. 

c. Random Errors vary from round to round in such an unpredic- 
table manner that no correction can be introduced by the fire control 
system. These random errors are classified as Round- to-Round errors and 
are the result of variations in ammunition rounds and characteristics of 
the weapon system. 

As can be seen by the error classification, the mean (p) of the dis- 
tribution of any round fired is determined by the fixed biases of the 
system. That is, the system bias offsets or shifts the position of the 
standard normal curve from the target center of mass to the center of 
impact, (fO The variable biases which are considered constant for the 
particular engagement and each burst within the engagement cause a dis- 
persion (o’ ) about the impact point. The random errors provide an index 
of the additional dispersion (o^) distance from the mean of the dis- 
tribution of any round fired in the burst. 

In the firing of an automatic cannon, it is compulsory to calculate 
the effects of several of the variable biases which have traditionally 
been listed only as a system variable bias. Brodkin 01 reported that 
if the calculated effects of these variable biases can be accurately 
measured that this provides a more representative impact point. This 
is extremely important in automatic cannon gunnery, since deflection of 
the muzzle and muzzle velocity changes can realistically be encountered 



1 ? 



in the burst of fire. As a result, if several variable biases are cal- 
culated they must be included as an additional fixed bias of the system 
and deleted as a source of variable bias. 

3. Miss Distance 

The point of impact at the target plane can be addressed by 
using the X and Y axes for reference . In most cases the gunner will 
have a line of sight to the target at the center of mass of the target. 
This point will be at coordinates 0,0 of the target plane. The offset 
from the target centerpoint when the projectile arrives at the target 
plane is due to fixed system biases and will be defined as the expected 
"Miss Distance." This expected miss distance is defined in both the 
azimuth (p-X) and elevation (p'Y) directions at the target. Synony- 
mously this miss distance is termed "X-error” and "Y-error" for azimuth 
and elevation respectively. (See Figure 1) 




Line of 
Sight 




Figure 1: Miss Distance Representation 
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This study attempts to mathematically model the miss distance 
for any round at the target plane by considering the system fixed biases, 
as well as accounting for any variable biases which can be accurately 
measured during the subject firing engagement. In order to compute 
X-error, the total displacement of the round from a given muzzle posi- 
tion must be determined. Dickinson W reports this lateral displace- 
ment from the line of fire is caused by crosswind, rotation of the earth 
and drift as addressed below. 



projectile along its intended route of flight. Most rounds which exhib- 
it drift are spin-stabilized and have high spin rates. Any drift de- 
flection is positive and arises from the clockwise spin of the projectile . 



as a result of Coriolis acceleration. However, the error is assumed for 
this study to fall into the category of Occasion- to-Occasion variable 
bias due to random characteristics. (<JX ? ,(JYg) 



flight will cause a deviation in azimuth from its intended trajectory. 
Normally this displacement is listed as a variable bias but it is felt 
that this must be addressed in the X-error computation for accuracy of 
trajectory under various Quasi-Combat Weather Conditions, since weather 
will be a definite factor influencing automatic gun firing. 



displacement of the round in the same direction at the target. This de- 
flection is caused primarily by tube or barrel bend, improper gun align- 
ment, and trunnion force effect of previously fired rounds in the 
burst. 



a. Drift (X.^) is caused by aerodynamic forces which act on the 



b. The rotation of the earth (Xg.Yg) causes a system fixed bias 



c. Any crosswind (X,^) force present during the projectile's 



d. Additionally, a deflection of the muzzle 
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In past analysis of hit probability, Jump (Xj) and Parallax (Xp) 
have been listed as a fixed bias of the gun system. However when moder- 
ately sophisticated fire control systems are employed, these errors are 

usually fully corrected according to Pfleger and Bibbero M- As such, 

2 

these fixed biases are considered a round dispersion (CT ) factor and 
will not be a factor in any X-error model for evaluation of advanced 
technology gun systems. 

The derivation of miss distance in the Y direction (Y-error) from 
the line of fire is caused by the following factors reported by Pfleger 
and Bibbero M. 

a. Deviation from projectile standard muzzle velocity (Y^) is 
the result of powder temperature variations, tube heat transfer and 
tube bend . 

b . Any deviation from the standard air density (Yp) will cause 
drag effects of the air mass through which the projectile passes toward 
the target. Accordingly, any percent change from standard air density 
due to meteorological effects will increase or decrease the range of the 
projectile . 

c. Deviation from standard air temperature (Y^) will cause an 
increase or decrease in range of projectile similar to air density 
changes . 

d. Any rangewind (Y,^) force will cause a deviation in elevation 
from the projectile's nominal trajectory. 

e. Range error fixed bias (Yp) is caused by inaccurate range 
estimation. Pfleger and Bibbero M report that for the laser range- 
finder this may be taken as zero, but when distance measurements are 
made by the human eye, the fixed bias cannot be stated with certainty. 
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Thus, range estimation should only be considered as a variable bias 
(CTY^) since the crew’ s estimation of the range will probably change in 
a random fashion for each engagement. 

f. Any elevation of the muzzle (Y^) deviating from the appro- 
priate superelevation setting for the specific range will increase or 
decrease the projectile range. This deviation is caused by gunner 
superelevation error, barrel bend, and trunnion force effect from previ- 
ously fired rounds in the burst of fire . 

As was previously noted in the X-error biases, there are certain 
biases which have traditionally been considered fixed but due to modern 
fire control systems, the assumption is made that they are fully correc- 
ted. Likewise, in computing Y-error at the target, the fixed biases of 
Jump (Yj) and Parallax (Yp) are considered to be zero and thus fully 
corrected . 

In this analysis the use of X-error and Y-error, miss distances, 
give the analyst a good indication of the point of impact for either a 
single shot or a burst of several rounds. Chapter HI, Section 3 in- 
cludes a discussion of the impact of barrel bend on the automatic gun 
and the resultant position of the muzzle at time of firing, at which 
time the approximate point of impact can be determined for each round. 

By knowing the mean (h>) of the standard normal curve determined by the 
fixed bias of the system and the dispersion about that point in the form 
of a variance (<7 ) , then integration of an assumed distribution over 
the target area is required to generate hit probabilities. 

4. Mathematical Model for Miss Distance 

The computation of miss distance requires a mathematical model 
which will express miss distance as a function of the sources of error 
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affecting the firing trajectory of the projectile. The variables which 
have an effect on miss distance at the target plane have previously- 
been enumerated in this chapter. Using these variables in conjunction 
with the appropriate firing tables for a particular projectile, the model 
can be derived to give a good approximation of miss distance at the 
particular range in question. 

This model is based on the following assumptions: 

a. Variable biases can be combined with the fixed system biases 
to refine impact point ( p- ) if the value can be measured directly and 
hence the effect determined. 

b. Jump (Xj,Yj) fixed bias equals zero due to modem fire con- 
trol systems and is treated only as an element of system dispersion (c ) . 

c. Parallax (&p,Yp) fixed bias equals zero due to modem fire 
control systems. 

d. Range Estimation (Y^) fixed bias equals zero and is addressed 
only as an element of the system dispersion 

e . The firing table data used is accurate for the particular 
round in question and will provide appropriate variable coefficients. 

f. The horizontal (X) and the vertical (y) errors are indepen- 
dent of each other. 

g. Fire control residual sensing error is accurately assessed. 

h. Deflection of the muzzle as a result of trunnion forces 
equals zero for any round fired. This is in fact not realistic but 
little research in this area is available for the study. 

The miss distance model for X-error (X„) is developed by con- 
sidering the variables which affect the lateral displacement of the 
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round from the line of fire . The variables previously listed in Chapter 
III are: 

a. Drift 

b. Cross-wind 

c. Muzzle lateral displacement resulting from inaccurate gunner 
laying, barrel bend, and effect of trunnion forces. 

The predicted point of impact of the projectile in meters at the 
target can be modeled using these variables in conjunction with firing 
table data for the round and weapon system as illustrated below. 

X T = C D + a l X W + ^ W 

where the coefficients a^ and a^ are given by: 

_ d Drift 

1 dm/s crosswind 

which is the change in projectile drift for a unit change in crosswind 
in meters per second from the firing table 

a _ d Lateral displacement 

2 " d Lateral muzzle position 

which is the change in lateral displacement of projectile for a unit 
change in lateral muzzle position expressed in mils. 

The coefficient, a. 0 , is computed by use of a simple geometric 
triangulation for direct fire weapons as shown in Figure 2. The width 
of the change can be computed by the following formula. 

Width (a 2 ) = MILS X 2gSSe (2) 

For example, at 1000 meters, a 1 mil (/) change in muzzle angle toward 
the target will result in a 1 meter lateral displacement of the round 
at the target plane . 
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Figure 2: Geometry of Muzzle Lateral Displacement 



The constant and variable derivations are: 

Cp = Drift of the round from firing table data at the desired 

range 



Xy. = Crosswind in meters per second 

(1) if cross wind is left to right = positive 

(2) if cross wind is right to left = negative 

= Muzzle lateral displacement in mils resulting from the 

summation of gunner error (X r ) , barrel bend (x^) , and trunnion force 

deflection (X^) . Thus: 
r 



*M = X G + S + X F 



( 3 ) 



This model for X-error is valid only if no fire control system 
wind sensors are employed. If sensors are employed then the model must 
be modified to account for the accuracy of the sensor since any correc- 
tion is a function of sensor and correction accuracy. The X-error will 
be the residual error remaining after sensing and corrective action to 
the system. 
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Past data on wind magnitudes indicates that an average of approx- 
imately 50 % wind sensing and correction is possible but this is again a 
function of the sensor accuracy. Thus for the system being evaluated, 
the sensing error must be calculated correctly and applied to Equation 2 
in the following manner. 



*T * °D + a h X w * V + W 

where is the approximate residual sensing error for crosswind. 

The miss distance model for Y-error (Y^,) is developed by con- 
sidering the variables which affect the vertical displacement of the 
round from the line of fire to center of target. These variables were 
previously listed in Chapter III ass 

a. Deviation from projectile standard muzzle velocity 

b. Deviation from standard air density 

c. Deviation from standard air temperature 

d. Range wind 

e. Muzzle vertical displacement resulting from inaccurate 
gunner laying and barrel bend. 

The model below can be used to estimate the Y-error in meters at the 
target plane by inserting these variables into the general equations 



Y T = - C 1 + Vv + a 2 Y D + a 3 Y A + Vw + a 5 Y M 



(5) 



where the constants and coefficients are defined as follows: 



-Cl - -(Superelevation^) * 



a _ d elevation 

1 d m/s standard muzzle velocity 
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which is the change in elevation for a unit change in meters per second 
from standard muzzle velocity as published for the projectile. 

_ d elevation 

2 ~ d percent standard air density 

which is the change in elevation for a 1 % change from standard air 
density as listed in the firing tables 

d elevation 

— ■■■ ■ . . 

3 d percent standard air temp 

which is the change in elevation for a 1 % change from standard air 
temperature as listed in the firing tables 

d elevation 

4 ~ d me ter/ sec. rangewind 

which is the change in elevation for a unit change of rangewind as 
listed in the firing tables 

d elevation 

3 , — 

5 d muzzle elevation 

which is the change in elevation for a 1 mil elevation of the muzzle 
per firing tables . 

Once the general model is established for the specified round and range, 
then Equation 5 is used to establish Y-error at the target plane by in- 
put to the following variables: 

Yy = the change in meters per second from standard muzzle 
velocity of the projectile specification data 

(1) if the muzzle velocity increases, then Yy is positive 

(2) if the muzzle velocity decreases, then Yy is negative 
Y^ = the percentage increase/decrease from the standard air 

density used in computing firing table data since firing tables are 
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computed using the International Civil Aviation Organization (ICAO) 
standard atmosphere . Standard density must be adjusted accordingly as 
altitude above sea level and geographic area are changed. The model for 
calculating the air density in gms/cc is defined in Humphreys jVJ as: 



1.2930 * 10 



-3 



B-0.378e 



1 + 2.039 * 10" 3 (t-32) 29-921 



( 6 ) 



where 

t = air temperature in degrees Farenheit 
B = total pressure in inches of mercury 
e = vapor pressure in mm of mercury 

By use of the tables in Smithsonian jj3j , a vapor pressure can be deter- 
mined for any temperature and relative humidity (RH) . Similarly, B 
(the total pressure in inches of mercury) can be determined at any 
standard atmospheric altitude from Humphreys [VJ . 

Y^ = the percentage increase/decrease from the standard air 
temperature . The ICAO standard air temperature at sea level (0 altitude) 
is 15°C and must be adjusted for altitude variation. 

Yy = rangewind in meters per second 

(1) if the direction of the wind is from the target toward 
the combat vehicle along the line of fire then Y^ is negative 

(2) if the direction of the wind is from the combat vehicle 
toward the target along the line of fire then Y u is positive . 

Y^ = the verticle muzzle displacement in mils (/i) resulting from 
gunner superelevation adjustment ( Y _ ) , barrel bend (Y_) , and trunnion 

U D 

force elevation (Y ) . Thus, 

F 

Y = Y + Y + Y (7) 

M G B F K(> 
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As with the mathematical model for computing miss distance in 
the X-direction, this model similarly is valid only when no sensors are 
used for sensing wind, air temperature and air density. When sensors 
are used, the residual error is a function of the sensing accuracy of 
the fire control system. For incorporation of sensors into the model, 
the percentage of each residual error source must be calculated. After 
calculation of these percentage values for the variables in question, 
the Y-error model becomes: 

Y t = -C t + a lY „ + a 2 (Y D * Sj,) + a,(Y A * S A ) + ^ 

V 1 * * V + Vm 

where S^, , and S T ^ are the percentage of residual sensing error for 

air density, air temperature, and rangewir.d respectively. 

5 . Round Dispersion 

3y using the mathematical models previously developed for both 
X and Y Miss Distance, the theoretic impact point at the target plane 
can be determined for any round fired either in a single shot or burst 
fire mode. The X and Y-error constitute an offset between the mean of 
the normal curve (center of impact) and the center of the target 
according to Pfleger and Bibber o M- It can be seen that the greater 
the expected miss distance due to fixed biases, the lower the hit 
probability. 

In Chapter III, Section A2 all classes of gunner/fire control 
system errors were defined. The magnitude of these errors - both fixed 
bias or miss distance, and variable errors - is affected by the sophis- 
tication of the fire control system used. 
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Since sophistication of the basic miss distance models has 



already been accomplished to account for fire control system product im- 
provement, the next area of concern is the total variation of the system 
about the target centerpoint of each round fired. Here again, it is felt 
that fire control system product improvement must be addressed . 



The ARPA HIMAG program is working with several baseline fire 



control system configurations to study the effects of component improve- 
ments and in defining minimum performance requirements. For purposes of 
this analysis it was felt that the model could be exercised by using the 
following two HIMAG baseline configurations as listed in General Motors 

Corporation GO. 



a. M60A1 Approximation: 

Slaved Ballistic Reticle Az,El 
Superelevation Only - No Sensors 
Manual Range Estimations 
Direct View Day Only OPTICS 

b. XM-1 Approximation: 

Independently Stabilized Gun Director Az,El 
Full Solution Gun Offsets 
Laser Rangefinder 
Sensor and Manual Estimates 
Direct View Optics (Day) , FIR Biocular 
Using the classes of error from Chapter III , Section 2 and the 



breakdown of these errors by class can be made . 

"Round- to-Round" errors (which are totally random since no 



theory of error components 




correction can be made to the system by the fire control system) are 
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caused by: 

a. Ammo Dispersion (C"X^,CTY^) 

b. Jump Dispersion (<XXj,<XYj) 

These random errors will be the same for any type fire control system 
and range. Round-to-round errors take on different values for each 
round fired in a burst and Cerrato D] considers them to be independent 
and normally distributed. The combined standard deviation in one dimen- 
sion for these random errors is calculated from the component error 
sources where 

0’*SH = (CrX A 2 +<7X /)* « 

and 

°-*hb = ( <ty a 2 +crY J 2 >* < 10 > 

"Burst- to-Burst" biases are random and by their very nature will 
take on a constant value for the entire burst of fire. These errors are 
caused by: 

aY L 

b. Visual Resolution (<T '/L^, (TY^) 

It is estimated by Cerrato bl that these errors will be independent and 
normally distributed and the combined standard deviations are of the 
form: 

<rx 3B = (crx 1 2 + c rx v 2 )^ (11) 

and 

CTY B3= ( crY l 2 + < ^ Y v 2 ) i <12) 



a. Gun Pointing (Laying) (CTX^, 
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"Occasion-to-Occasion" biases change so slowly that they can be 
considered constant for the period of engagement against a stationary 
target. Falling into this category are: 

,CTY b 

b. Fire Control (<XX F ,<XY F ) 

c. Ranging (CTY R ) 

d. Cant (CTX c ,(TY c ) 

e. Angle of Site (cTY<.) 

f. Earth Rate (cJXg.CT Yg) 

g. Tube Droop (<T Y^) 

h. Zeroing (1200 M) (jX^C Y z ) 

For approximation purposes these errors can be considered independent 
and normally distributed, and the combined standard deviations are of 
the form: 

<T x Q0 = (c rx B 2 +<tx f 2 + crx c 2 + crXg 2 + c rx z 2 )* ( 13 ) 

and 

<^00 = (°' Y B 2 +<rY F 2 +CTY R 2 tCTY C 2 +<JY S 2 + ,,,, 

9 9 9 1 . 

+<7Y D +<XY Z 

These errors are not presented as an all inclusive list for 
every weapon system in use . They represent the authors ' best list for 
use in approximating the fire control systems being presented in this 
study. Additionally, the computations performed in computing the in- 
dividual standard deviation values are based on guidelines established 
in Pfleger and Bibbero [VJ and are approximations only. Accordingly, 
more exact values should be used and additional errors accounted for 
when available for use . 



a. Boresight Loss (dX^ 
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After categorizing the various errors for analysis, an approxi- 
mation of the system error magnitudes for total dispersion about the miss 
distance point can be computed. Tables I and II were generated by 
combining assumed firing table data with suggested error magnitude 
computations in Pfleger and Bibbero 03 for both the ARPA APFSDS and 
HE rounds . 

6 . Model Application 

Using assumed firing table data for both the ARPA APFSDS and HE 
rounds listed in Appendix A, the miss distance model can be applied to 
both the M60A1 and XM-1 fire control system configurations. By gener- 
ating the constant and coefficients for the model, these equations can 
then be used in the simulation of miss distance presented in Chapter 
IY. 

To compute X-error, the mathematical models developed in 
Chapter III, Section A4 are used. (See Equation 1 for M60 FCS con- 
figuration and Equation 4 for XM-1 FCS approximation) . Inserting 
appropriate firing table data at the desired range into the model re- 
sults in Miss Distance formulas for use in simulation. These formulas 
are summarized in Table HI and IY. 

The mathematical models previously developed in Chapter HI, 
Section A4, are used for computing Y-error at the target plane. (See 
Equation 5 for M60 FCS configuration and Equation 8 for XM-1 FCS approxi- 
mation) . Tables V and VI were generated for the ARPA APFSDS and HE pro- 
jectiles by inserting assumed firing table data into the basic mathe- 
matical models. Additionally it is assumed that residual sensing errors 
are 1 %, 2 % and 3 % respectively for density, air temperature, and 
range wind . 
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TABLE I 

FIRE CONTROL SYSTEM ERROR MAGNITUDES 
APFSDS ROUND 
QUASI -COMBAT CONDITIONS 





M60A1 APPROXIMATION 


XM-1 APPROXIMATION 


1500 M 


2500 M 


1500 M 


2500 M 


Error Source 


Azimuth 


Elev 


Azimuth 


Elev 


Azimuth 


Elev 


Azimuth 


Elev 


Round-to- 

Round 




<rx A!_«iL 


.320 


.280 


.320 


.280 


.320 


.280 


.320 


.280 




.220 


.300 


.220 


.300 


.220 


• 300 


.220 


.300 


Burst- to- 
Burst 






.254 


.254 


.153 


•153 


.020 


.070 


.020 


.070 


<y\, ^ 


.060 


.060 


.060 


.060 


.060 


.060 


.060 


.060 


Occasion- to- 
Occasion 






.080 


.080 


.080 


.080 


.080 


.080 


.080 


.080 




. 0?0 


.070 


.070 


.070 


.070 


.070 


.070 


.070 




- 


.624 


- 


1.120 


- 


.006 


- 


.006 


n c (n c 


.252 


.001 


.436 


.001 


- 


- 


_ 


_ 


^s 


- 


.011 


_ 


.019 


_ 


.011 


_ 


.019 


«r<™E 


.075 


.075 


.125 


.125 


.075 


.075 


.125 


.125 


a h’ a \ 


.212 


.220 


.294 


.242 


.170 


.185 


.198 


.210 



NOTESs 



(l) All error values represent one standard deviation (l ) of the error 
source measured in mils 



(2) Xl = * 



1000 

R 



y l = 



3 



where a.^ and a^ = •§■ target size in each dimension (mils) 

R = Target Range (meters) 

(3) Gun laying errors in XM-1 FCS configuration reflect testing results 
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TABLE II 

FIRE CONTROL SYSTEM ERROR MAGNITUDES 
HE ROUND 

QUASI-COMBAT CONDITIONS 





M 60 A 1 APPROXIMATION 


XM -1 APPROXIMATION 


1500 M 


2500 M 


1500 M 


2500 M 


Error Source 


Azimuth 


Elev 


Azimuth 


Elev 


Azimuth 


Elev 


Azimuth 


Elev 


Round-to- 

Round 




^A’^A 


.320 


.280 


.320 


.280 


.320 


.280 


.320 


.280 




.220 


.300 


.220 


.300 


.220 


.300 


.220 


.300 


Burst- to - 
Burst 




c \ 


.254 


.254 


.153 


•153 


.020 


.070 


.020 


.070 


ffX I’ aI 2 


.060 


.060 


.060 


.060 


.060 


.060 


.060 


.060 


Occasion- to- 
Occasion 






.080 


.080 


.080 


.080 


.080 


.080 


.080 


.080 


m f 


.070 


.070 


.070 


.070 


.070 


.070 


.070 


.070 


^r 


- 


1.684 


- 


3-792 


- 


.016 


- 


.022 




.591 


.001 


1.144 


.001 


- 


- 




- 


- * Y S__ _ 


- 


.025 


- 


.049 


- 


.025 


- 


.049 


*V «E 


.075 


.075 


.125 


.125 


.075 


.075 


.125 


.125 


<rX Z' aT Z 


.341 


.220 


.605 


.242 


.170 


.185 


.198 


.211 



NOTES: 

(l) All error values represent one standard deviation (l ) of the error 
source measured in mils 



(2) 

where a^ and = -§- target size in each dimension (mils) 

R = Target Range (meters) 

( 3 ) Gun laying errors in XM -1 FCS configuration reflect testing results 
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TABLE III 

MISS DISTANCE X-ERROR FORMULAS 
APFSDS ROUND UNDER QUASI -COMBAT CONDITIONS 

M60A1 FCS APPROXIMATION 
1500 Meters 

X T = 0 + .0358X W + 1.5^ 

2500 Meters 

X T = 0 + .0983X w + 2 . 5 \ 

XM-1 FCS APPROXIMATION 
1500 Meters 

x T = 0 + .0358(x w * .03) + 1.5^ 

2500 Meters 

X T = 0 + .0983(X W * .03 ) + 2.5X m 
NOTE: Residual Sensing Error for Crosswind Assumed = 3 % 
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TABLE IV 

MISS DISTANCE X-ERROR FORMULA 
HE ROUND UNDER QUASI-COMBAT CONDITIONS 

M60A1 FCS APPROXIMATION 
1500 Meters 

Xj = 0 + •197X W + 1.5^ 

2500 Meters 

X T = O + .614^ + 2.5X m 

XM-1 FCS APPROXIMATION 
1500 Meters 

X T = 0 + .197(3^ * .03) + 1 . 5 \ 

2500 Meters 

X T = 0 + .614(X W * .03) + 2 . 5 \ 

NOTE; Residual Sensing Error for Crosswind Assumed = 3 % 
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TABLE V 

MISS DISTANCE Y -ERROR FORMULA 
APFSDS ROUND UNDER QUASI-COMBAT CONDITIONS 



M60A1 FCS APPROXIMATION 



1500 Meters 



Y t = -4.179 + .006Y V + (-.0025Y d ) + (-.001JY a ) + 
•OOOY m + 1.473Y m 

2500 Meters 

Y ? = -12.007 + .0l44Y y + (-.011Y d ) + (-.003Y a ) + 
.ooi5Y w + 2.454 y m 

XM-1 FCS APPROXIMATION 



1500 Meters 



y t = -4.179 + .006Y y + -.0025 (Y d * .01) + 

-.0015(Y a * .02) + ,000(Y W * .03) + 1.473Y m 

2300 Meters 



Y 



T 



-12.007 + .0l44Y y + 
-.003 (Y a * .02) + 



-.oii(y d * .01) + 

.0015(Y W * .03) + 2.454Y m 
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TABLE VI 

MISS DISTANCE Y -ERROR FORMULA 
HE ROUND UNDER QUASI -COMBAT CONDITIONS 



M 6 OAI FCS APPROXIMATION 
1500 Meters 



Y t = - 9.7766 + .0198Y y + (-.0218Y D ) + (-.008Y A ) + 
•0033Y W + i.4738Y M 

2500 Meters 



Y t = -31.5553 + -0633Y v + (-.1298Y D ) + (-.0353Y A ) + 
.0153Y V + 2.4553Y m 

XM-1 FCS APPROXIMATION 



1500 Meters 



Y t = - 9.7766 + .0198Y y + -.0218 (Y d * .01) + 

-.008(Y a * . 02 ) + .0033(Y W * .03) + 1.4?38Y m 

2500 Meters 

Y T = -31.5553+ -0633Y V + -.1298(Y d * .01) + 

-.0353(Y a * .02) + .0153(Y W * . 03 ) + 2.4553Y m 
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7. Summary 



This section has presented a basic Miss Distance model that can 
be used to determine the impact point of each round fired at the target 
plane, either in single shot or burst mode of fire. Hit probabilities 
(P„) have not been computed because miss distance is the basis for any 

n 

later generation of P„. Although the model as presented is simplistic 

n 

in its formulation, since much of the laborious computation is aleviated 
by using published firing table data, the data calculations for insertion 
into the model variables are not so readily available. This is espe- 
cially true when attempting to numerically define tube (barrel) bend, 
effect on the muzzle due to trunnion forces, and nonstandard muzzle 
velocity changes due to nonstandard propellant temperature and barrel 
bend. The examination of all of these input variables is beyond the 
scope of this study, however, the remaining portion of this chapter will 
be dedicated to an in depth study of barrel bend and its resultant im- 
pact demonstrated in Chapter IY, The Simulation Model. 
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B. BARREL BEND 



1 . Introduction 

It has long been observed that if a strip of metal is heated and 
then allowed to cool in such a manner that one side cools more quickly 
than the other, the strip will bend toward the cooler side. A more 
accurate description of this process would be to say that the cooler 
side returns to its original temperature and shape more quickly. 

This process also occurs in gun barrels. When the gun is fired, the 
barrel is heated. The barrel is exposed to the environment as it cools. 
If sides of the barrel are subjected to varying environmental conditions, 
e.g. wind or rain, the barrel will cool unevenly. This uneven cooling 
causes a temperature differential across the barrel. The magnitude of 
this temperature differential increases over time to a maximum value 
and then diminishes as ambient air temperature is approached, as illus- 
trated in Figure 3* 



Temperature 

T 

am 



Time 

T q = Temperature of hotter side 

T^ = Temperature of cooler side 

T = Ambient air temperature 
am * 

Figure Time -Temperature Curves 
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To visualize the establishment of this differential more clearly, 
imagine that the temperature of all points of a gun barrel cross section 
has been raised to a particular value. Suppose further that a wind is 
blowing at point A in the figure below. 




Figure 4: Gun Barrel Cross Section 



Point A will receive the full benefit of the cooling effect of 
the wind and cool much more quickly than Point C. Moving from Point A 
toward Points B or D the temperature will rise until, in the vicinity 
of those points, it equals that of Point C. The temperature would then 
remain relatively constant until Point C is reached. A temperature 
curve for the surface would appear as shown in Figure 5* 

Rain will cause a similar differential. In this case the hottest 
point would likely be near Point D and the coolest near Point B. The 
sun can also affect a temperature differential by adding temperature to 
the top of the barrel. 

If there is a temperature differential along a gun barrel, the 
barrel will bend in the direction of the cooler side. If a round is 
fired through a bent gun barrel there will be an effect on the intended 
trajectory of the round. In order to approach the problem of 
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Temperature 




Surface Points 

T q = Temperature at Point C 
T q -AT = Temperature at Point A 

Figure 5: Temperature Curve 



determining the effects on accuracy, it is necessary to determine how 
hot the outer surface of the barrel gets, how much of a differential can 
occur across the barrel, and how much bend will result. 



2. Approximation of Gun Barrel Surface Temperature (T ) 

The highest temperature on the outside surface of a gun barrel 
(T ) is a complex function of many factors including barrel thickness, 
material, type of propellant, number of rounds fired, rate of fire, etc. 
T q is also unique for each different type of gun and will vary slightly 
between similar guns. Heat and temperature curves for many guns have 
been developed. Since the 75ram cannon has been selected for this 
analysis, T q must be approximated for this particular barrel. 

Figure 6 shows temperature curves developed for temperature as 
a function of rounds fired at maximum rate of the ARES 75"® cannon. fro] 
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Figure 6: APES Temperature Curves 
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The area of concern is near the muzzle end of the barrel and can be 
assumed to be representative of the temperature along the forward end of 
the barrel. The outer barrel temperature is the quantity of interest 
for the analysis. The curve in Figure 6 was developed by a computer 
program and the exact mathematical basis could not be provided. The 
shape of the curve indicates that a linear approximation should be very 
satisfactory for estimating T q in the region of interest. Linear re- 
gression on data points of the curve results in the following approx- 
imation for T . 

o 

T q = 64.55 + 19. 8X (15) 



where X = number of rounds fired at maximum rate . 



Figure 7 shows the data from the Figure 6 curve in the first 



two columns. T 

0 


estimated 


from (15) is shown in the 


third column. In 


the last column 


is T q estimated by substituting the 


ambient air temper 


ature (T ) for 
' am' 


the value 


of 64 .55 in ( 15 ). 




Rounds 




Regression 


Regression 


Fired 


Figure 6 


(with 64.55) 


(with T ) 
-» am' 


0 


75 


64.55 


75 


5 


125 


I 63 .OO 


174 


10 


225 


262.00 


273 


15 


360 


362.00 


372 


20 


485 


460.00 


471 


25 


595 


559.00 


57 0 


30 


690 


658.00 


669 


35 


780 


757.00 


768 


40 


860 


856.00 


867 


45 


940 


956.00 


966 


50 


1010 


IO 55 .OO 


IO 65 



Figure 7: Approximation of T q 
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The original regression and the regression with the substitution 

of T are adequate approximations, especially in the range of 10 to 45 
am 

rounds. This approximation is used to estimate the temperature differ- 
ential in the development of subsequent relationships. Trials with 
differing ambient temperatures show that while T q varies, as would be 
expected , the temperature differential remains the same . Also the use 
of the approximation with T am included simplifies some later calcu- 
lations. The coefficient of determination (R ) of the regression in- 
cluding T is greater than .99. Therefore, for the balance of the com- 
am 

putations made , T q will be approximated by the following: 



T = ambient air temperature 
am r 

X = number of rounds fired at maximum rate 

3. Approximating the Temperature Differential 

A formula for the effect of wind or rain will take the general 



T = T + 19. 8X 



o am 



( 16 ) 



where 



form: 



AT = C(T - T ) 0<C < 1 

o am 



( 17 ) 



where 



AT = the temperature differential 



C = cooling factor 



T q = temperature of the hottest part of barrel surface 

T = ambient air temperature 
am 
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In other words, the temperature differential that will develop 
across the barrel will be some fraction of the difference between the 
hottest point of the barrel and the ambient air temperature. The cool- 
ing factor, C, will be the value attainable by some function of the 
intensity of wind or rain, time and other factors. The search for 
functions representing the cooling effect of wind or rain was most un- 
successful. Neither functions nor sufficient data to permit derivation 
could be found. The search was then redirected to find the maximum ex- 
pected value for C. A value of G could not be determined for rain. 
However, determination of the maximum expected value of G is possible 
for wind cooling. 

By Newton's Law of Cooling: 

Q = h A AT (18) 



where 

Q, = heat transfer rate 

h = heat transfer coefficient 

AT = temperature differential between the 
surface and ambient air temperature 

We can reasonably assume the heat transfer rate (Q) on the sur- 
face of a gun barrel to be constant. The definition of the temperature 
differential is: 



AT = (T - T ) 
am 7 



(19) 



We can combine (18) and (19) under the above assumption to give: 



(T - T ) h = 
x am' 



_ Q. 
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Q 

A 



( 20 ) 




Figure 8 s Gun Barrel Gross Section 



Since the temperature differential "between and T am in 
Figure 8 is different from the temperature differential between T q and 
Tam and since both quantities equal a constant value, then the heat 
transfer coefficient (h) must vary and: 



(T. - T ) h = -v- 
1 am' max A 

(T - T ) h . = - 7 - 

v 0 am' min A 

From (2l) and (22): 

(T, - T ) h = (T - T ) h . 

1 am' max o am' min 

Since, by definition: 



T, = T - AT 
1 o 



( 23 ) can be reduced to: 

h 



T _ T _ max 

o am h . 

min 



(T - T - £T) 
' o am ' 



( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 
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Solving (25) for AT yields: 



h . 

AT = (1 - r~-) (T - T ) (26) 

' h ' v o am' v ' 

max 

The ratio of h . and h can be determined from the work of 
mm max 

Giedt . [u] With a Reynolds number of 70,800 and 4.25 inch diameter gun 
barrel [ 12 ] the value of the ratio is found to be 100/250 or 0.4. Thus, 
the wind cooling factor for the ARES cannon is 0.6. This corresponds 
roughly to a 20 mph wind. Formula (26) for the temperature differential 
may now be reduced to: 

AT = 0.6(T - T ) (2?) 

w ' o am' ' 

where 



AT = the maximum expected temperature differential 
across the barrel of the ARES cannon caused 
by wind cooling 

T q = temperature of hottest point of barrel surface 

T = ambient air temperature 
am p 

With the above expression for the maximum expected temperature 
differential, it is possible to calculate the magnitude of bend that can 
occur in the ARES cannon as a result of wind cooling. This requires a 
formula for bend which is developed in the following section. Factually, 
this is the only bend that can be calculated with some confidence. It 
is desirable, however, to have some feel for the range of bend that can 
occur because of varying conditions of rain, as well as wind. There- 
fore, a parametric analysis was performed with assumed values for other 
conditions. It must be emphasized that the value of 0.6 for wind cool- 
ing is the only value that has been derived. This value will be used in 



the study as the cooling factor for "hard wind." Other values are based 
solely on the subjective judgement of the author and others consulted. 



For the purpose of the parametric study the following cooling factor 
values will be assumed for the weather conditions indicated t 



Haxd wind - 0.6 
Light wind - 0.2 
Hard rain - 0.8 
Light rain - 0.6 



(approximately 20 mph) 

(approximately 5 mph) 

(approximately 0.5 inches per hour) 
(approximately 0.05 inches per hour) 



4. Derivation of the Barrel Bend Formula 



From the definition of curvature: 



K = 



dQ 

dx 



where 



K = curvature 
9 = local deviation angle 
x = desired axis 
it can be seen that 




(28) 



(29) 



The task is to find an expression for K. Hooke’s Law for isotropic 
materials 13 states: 





(30) 



where 



£ = strain in direction of length 



R. S. Newton and M. D. Kelleher of the Department of Mechanical 
Engineering, Naval Postgraduate School provided assistance here and in 
other parts of the derivation. 
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CT x = normal stress on a barrel cross section 
E = modulus of elasticity 
CC = coefficient of thermal expansion 
T = temperature of a point 
T a = average temperature 
Multiplying (30) by y da gives 

JV X y da = jV x y da + <*£(T - T a ) y da ( 31 ) 

The bending moment does not change as a result of temperature strains 
and may be considered zero for this purpose. Therefore, 

^-jV x yda=0 ( 32 ) 

and ( 31 ) may be reduced to: 

jV x yda= q^(T - T a ) y da (33) 

Since by Popov [13] 

i x = K y (3^) 

(33) becomes 

K Jy 2 da = tfJ(T - T a ) y da (35) 

By definition: 

y 2 da = I = (r ^ - r.^) ( 36 ) 

j z 4 v o 1 ' 

where 

r = radius of outer surface 
o 

r. = radius of inner surface 
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Substitution of ( 36 ) into (35) yields: 



K = C (T - T a ) y da ( 37 ) 

z J 

It is necessary to make certain statements and assumptions before 
proceeding with the derivation. If a gun barrel cross section is ex- 
periencing uneven cooling, a temperature differential will develop. As 
illustrated in Figure 9* the surface temperature at 0 = 0 will be equal 
to T q -AT. The temperature will rise as the angle 0 increases and 
reach a maximum value near 0 = The temperature will remain con- 

stant until 0 = TC. Assume the point of equal temperature occurs where 











Figure 9: 


Cross Section Temperature 




Then the surface 


temperature 


(T„) can be expressed by the 


following: 


T = T 
s 0 


- AT cos 0 


0 -0 < 


(38) 


T = T 
s 0 






(39) 
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The temperature distribution in a hollow cylinder is known by Krieth 
to be : 




T. - T 

T = T. - — In (— ) (40) 

i r v r. ' v ' 

In <£) 

1 

Finally, it can be shown that relatively little bending occurs 
at the breech end of the gun barrel. This fact is due to the increased 
wall thickness and because of the various supports that are found in this 
portion of the barrel. The amount of bend in this portion is assumed to 
be insignificant. Therefore, the portion of the barrel to be considered 
is the portion which extends forward of the supported section. This 
portion will be referred to as the unsupported length. It is also noted 
that in the unsupported length of the barrel, the inner and outer radii 
are generally constant for most guns, including the ARES cannon. If the 
radii are not constant in this region, then an expression must be in- 
serted in the relationships which define the radius as a function of 
length. Inner and outer radii will be considered constant through the 
unsupported length of the barrel for the remainder of the derivation. 

From (37) let: 

P = j* (T - T a ) y da (41) 



in cylindrical coordinates 
/'vTr r 



P = 



0 Jr. 

w l 



(T - T ) y r dr d0 

3 . 



(42) 



Substituting the definition y - r cos and (40) into (42) yields: 
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p = 



1° [ T i- (T i- 



v^r p r 
*0 



T ) 
s 7 



“ft) 



“ft 



- T^j cos 0 dr d0 



(43) 



Integration with respect to r gives: 
rtf' 



P = 



(k 



T ) (r ' 
a 7 ' o 



r. J ) - (T. - T ) 
1 J x 1 s' 



In 



r. 

^ i 



Tft-J Ht 



1 I r o 



9 \ r . 



cos 



0 d0 



(44) 



Let: 



T' = (T. - T ) (r 3 - r. 3 ) 
Y i a' v o i 7 



(45) 






In I 



r. 

l 



T&-) *1? 



1 r 0 



r. 

i 






Substituting (45) and (46) into (44) gives 

J dt' pit 

'X cos 0 d0 -X\ (T. - T ) cos 0 d 0 

0 J 0 1 s 



(4?) 



From the definitions of T g given by (38) and (39) > (47) may be expanded 



to: 



P = f J 'X cos 0 d0 - X C T. cos 0 d0 + A f (T -^T cos 0) 
JO JO 1 J 0 ° 



cos 0 d0 + A f T cos 0 d0 
¥ ° 



(48) 
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Integration with respect to 0 gives 



P = 



A at 
— u 



m 



(49) can now be replaced in ( 37 ) to give 

T, CL l A a t \ 

l -5 ) 

which by ( 36 ) can be further reduced to 



T , A of At 
K = zr 



4 

r - r. 
o 1 



( 51 ) can now be substituted into (29) to give: 

rx A a At 



e = 



dx 



0 r 



(50) 



(51) 



(52) 



The result of integration of ( 32 ) will be in radians . Since it 
is desired to have the result in mils, the right side of ( 52 ) is multi- 
plied by the conversion factor of 1018.6 mils per radian. Since the 
terms in ( 52 ) axe constants, including the mil conversion factor, let 
a bending constant (^) be defined where: 



6 = (ioi8.6)(X)(a) 

4 4 



Then (52) may be reduced to 



0 = 




£AT dx 



(53) 



(54) 



Since the terms to be integrated are constant along the unsupported 
length of the barrel, the result of integration will give the amount of 
bend at the muzzle in mils and will be: 



9 = L p AT 



(55) 
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where : 

L = length of unsupported barrel in millimeters 
P = bending constant defined by (53) 

AT = temperature differential defined by (26) 

5* Calculating Bend for the ARES Cannon 

Relationships derived in the foregoing sections permit calcula- 
tion of the maximum expected bend for the ARES gun as a result of wind 
cooling. For the ARES cannon [l2] : 
r Q = 54- mm 
r i = 37-5 mm 

L = 1700 mm 
OL = 6 x 10" 6 

Inserting the values above into (46) , (53) and (55) gives the final 
result: 



e = 0.0327 AT (56) 

(27) shows that for wind cooling: 

AT = 0.6(T - T ) (57) 

x o am' 

By (16) , T q may be approximated by: 

T = T + 19.8X (58) 

o am 

Substituting (58) into (57) and solving for AT, and then substituting 
the result into (56) yields 

0 = 0.39X (59) 

where X is the number of rounds fired at maximum rate. 
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This leads to the conclusion that the maximum expected bend at 
the muzzle of the ARES cannon as a result of wind cooling can be approx- 
imated by (59) . Recall that (58) is a good approximation especially 
in the range of 10 to 45 rounds. A graphic representation of (59) is 
given in Figure 10 . 



Bend 

(mils) 




Figure 10: Wind Generated Bend (ARES Gun) 



As can be seen from Figure 10, the amount of bend that can be 
caused by wind cooling can be very significant. These results are 
markedly less severe than has been estimated by other researchers. Rain 
can cause even more severe bend. 

It must be recognized that it is extremely unlikely that a tank 
crew would fire 45 rounds at maximum rate. However, given the capabil- 
ity, it is not unlikely that 10 rounds would be fired. Several shorter 
bursts fired in rapid succession will generate nearly the same amount 
of temperature and the resultant bend. An additional question must now 
be addressed. Given that the barrel is bent, how much deflection to the 
projectile will result? 
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6 . Straightening Forces 



If a round is fired through a bent barrel, the projectile may 
exert a straightening force on the barrel as it travels down the bore . 

If there is a degree of straightening, the amount of actual deflection 
will be less than the full amount of bend. Any amount of straightening 
that may occur must be paid for by an accompanying loss in muzzle veloc- 
ity. A loss in muzzle velocity will cause a range error. The round will 
not travel as far as intended and will pass through the target plane at 
a point lower than planned. There is no easy way to determine the net 
deflection and any change in muzzle velocity. 

The results of the parametric analysis conducted in Chapter TV 
indicate that the amount of bend will seldom exceed one degree. This 
amount of bend is very significant as far as accuracy is concerned, but 
relatively insignificant with respect to the forces involved. It is 
assumed that for the small angles encountered, the amount of straignten- 
ing will be negligible. The full amount of bend will then be considered 
as deflection of the round from its intended trajectory and there will 
be no resulting loss of muzzle velocity . 

7 . Time Dependency 

Figure 3 indicates the magnitude of bend varied over time. The 
angle will increase to a maximum magnitude and then diminish as each 
side of the barrel approaches the ambient air temperature. The maximum 
expected bend from wind cooling will occur approximately 15 minutes 
after completion of firing. A hard rain will cause the maximum bend to 
occur in a much shorter time . This time cannot be accurately calculated 
due to the lack of adequate data. It is assumed in this study that the 
maximum magnitude of rain induced bend will occur approximately five 
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minutes after completion of firing. The time required, for the maximum 
value of bend to occur will vary between types of guns and with varying 
amounts of wind or rain. 

8 . Avoiding Barrel Bend 

It is a fact that gun barrels bend. The foregoing sections have 
demonstrated that the amount of bend can be very significant and can 
cause a sizable miss distance at the target plane. This effect can be 
avoided by two general methods. The amount of bend occuring at any 
point in time can be compensated for or barrel bend can be eliminated 
altogether . 

Methods of compensating for barrel bend include: 

a. Temperature sensing of points of the barrel to determine 
temperature differential, calculating the amount of bend, and adjusting 
the lay of the gun accordingly . 

b. Direct measurement of bend through optic or other sensing 
devices and compensating as above . 

There are several ways of avoiding barrel bend altogether. 

These basically involve the prevention of uneven cooling. Two methods 
are: 

a. Installation of a cooling system which insures even cooling 
of all parts of the barrel. 

b. Shield the barrel from the causes of uneven cooling (i.e. the 
environment) with a protective shroud or coating. This alternative is 
reportedly being developed by certain allied nations . 

9 . Summary 

Gun barrels are heated during firing. Rapid fire guns in partic- 
ular can generate rather high barrel temperature, although high barrel 
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surface temperatures are certainly common in guns with slower rates of 
fire when a large number of rounds are fired. The gun barrel is ex- 
posed to varying environmental effects as it cools after completion of 
firing. Wind or rain on the surface of the barrel during the cooling 
period will cause the barrel to cool unevenly, resulting in a bend of 
the barrel. This occurrence can significantly affect the trajectory of 
rounds fired while the barrel is bent causing a sizable miss distance 
at the plane of the target. The effects of barrel bend can be compen- 
sated for or barrel bend can be avoided altogether. Weapons under 
development should be required to incorporate methods of avoiding barrel 
bend or its effects and modifications to existing weapons should be 
investigated . 
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IV. THE SIMULATION MODEL 



A. MODEL DESCRIPTION AND ASSUMPTIONS 

To employ the miss distance and barrel bend models under Quasi- 
Combat conditions, a computerized parametric analysis was conducted. It 
was performed in order to evaluate the automatic gun concept under 
various weather conditions, fire control system configurations and pro- 
jectile types. When required, the analysis used specific characteristics 
of the ARES 75rom automatic gun. 

To evaluate the effects of weather on gun performance (and in par- 
ticular the effects of barrel bend) a weather case matrix was developed 
containing variable inputs of air temperature, direct left to right 
crosswind, rainfall, vapor pressure, and cooling factors which were 
developed in Chapter III, Section B. The weather condition for any 
trial is considered to occur at 1000 feet above sea level. Air temper- 
ature varies from the standard of 60° Fahrenheit to either 100°F or 0°F. 
Rangewind has not been simulated but direct left to right crosswind 
effects can be determined by using a light wind of approximately 7 mph 
and a hard wind of approximately 20 mph. The effects of rain are simu- 
lated by considering a light rain condition of 0.05 inches per hour and 
a hard rain condition of 0.5 inches per hour. Vapor pressure values 
(in inches of mercury) are a function of the air temperature and relative 
humidity and are used in the model to compute air density. For the "no 
rain" condition, the relative humidity was assumed to be 80%. If a rain 
condition existed, then the relative humidity was set at 100%. By 
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referring to vapor pressure tables in Smithsonian Institute isj, spe- 
cific values were then generated for vapor pressure considering the air 
temperature and relative humidity. 

Each test weather case was conducted in conjunction with a selected 
projectile (HE or APFSDS) and range to engage the target. Once all 
selected weather cases had been simulated for a specific test projectile 
and range, then all weather cases were repeated with the next test pro- 
jectile and range . 

After generating the basic situation by establishing a weather case, 
the test projectile and range of the target, the only element remaining 
in the model scenario is the trial itself. The basis for the structure 
of a trial is an attempt to structure a somewhat realistic combat 
scenario for a rapid fire gun mounted on a tank. 

The trial represents perhaps an hour in combat for a tank. It is 
envisioned that during this period a tank would have several short 
engagements . During each engagement the tank would fire several short 
bursts of either three or five rounds at the target. Between engage- 
ments the tank could remain in position or move to a new location. For 
this study three engagements per trial and four bursts per engagement 
were used. (See Figure 11) 

The time between engagements is established as the time required for 
the maximum expected bend to occur. This time will vary according to 
the environmental conditions being experienced by the barrel. Maximum 
bend can occur in as little as five minutes under a hard rain condition 
or longer for lessor weather conditions. The length of any engagement 
is less than one minute . An engagement is thus considered a sequence 
of four short bursts fired in rapid succession at maximum rate of fire 
with minimal time between bursts. 
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ENGAGEMENT 



BURST 



ROUNDS/BURST 







3 or 5 
3 or 5 
3 or 5 
3 or 5 
3 or 5 
3 or 5 
3 or 5 
3 or 5 

3 or 5 
3 or 5 
3 or 5 

3 or 5 



Figure 11: Trial Representation 

The gun barrel surface heat generated in an engagement will nearly 
approach that generated by firing a single burst of the same number of 
rounds. The program structure calls for the gun barrel surface to be 
at ambient air temperature for initialization of each trial. Therefore, 
there is no barrel bend for the first engagement. The number of rounds 
fired during the first engagement is used to calculate the resulting 
gun barrel surface temperature . Bend at the muzzle in mils is then 
calculated depending on the weather case selected. Since wind is 
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assumed to be blowing left to right and perpendicular to the barrel, 
bend in the X axis is always to the left. Rain is assumed to be falling 
on top of the barrel. Rain induced bend is calculated in exactly the 
same manner with only the cooling factor changed. Thus, bend in the Y 
direction is always upward. 

For the third engagement of a trial one final factor must be con- 
sidered . The number of rounds in the second engagement generate barrel 
surface temperature which is added to the temperature remaining from the 
first engagement. Thus the bend generated between the second and third 
engagements will be more severe . After consideration of the time span 
and the foregoing assumptions, it was determined that a standard cumu- 
lative heat value of 0.5 should be used. Therefore, the barrel surface 
temperature at the end of the second engagement is based on the number 
of rounds in the second engagement plus 0.5 times the number of rounds 
in the first engagement . 

B. MODEL RESULTS 

Output tables from the computer trials are given in Appendix B. In 
order to simulate the two systems, the standard deviations developed in 
Chapter III, Section 5 were used. The basic differences in the two 
system approximations were due primarily to the higher sophistication of 
the XM-1 fire control employing sensor devices and the laser rangefinder. 

Table VII is a sample of the computer output to demonstrate the 
model’s capabilities. One trial is presented and the applicable general 
trial situation is listed in the page heading. In this case, the trial 
situation occurs using the ARPA HE round and engaging a target at a 
range of 2500 meters under weather conditions consisting of an air 
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TABLE VII 

COMPUTER OUTPUT SAMPLE 
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temperature of 0°F with a hard crosswind of approximately 20 mph and no 
rain. For reader clarity, some of the column headings presented for the 
trial results are defined below. 

1. FIXED BIAS (MILS) is the deviation in mils at the muzzle of the 
fixed biases of the weapon system for any round fired. 

2. BEND is the amount of barrel bend on the muzzle expressed in mils 
under the scenario conditions. 

3. TOTAL is the total fixed bias in mils at the muzzle under the 
scenario conditions, including barrel bend. 

4. RANDOM ERROR (MILS) is the effect of Round- to-R ound , Burst- to-Burst, 
and Occasion- to-Occasion errors on the projectile expressed in mils at 
the muzzle . These error values were generated using the error classi- 
fications and computation of standard deviations as presented in Chapter 
III and then applying the error values by use of a standard normal ran- 
dom number generator routine. After completion of the computer runs, a 
potential problem with the random number routine was discovered. There- 
fore, this routine should be completely verified before using the model 
in a production mode . 

5. AZIMUTH refers to the X axis of the target plane. 

6. ELEVATION refers to the Y axis of the target plane. 

7. TOTAL ERROR (METERS) is the miss distance in meters of the projec- 
tile at the target due to the contribution of all the fixed biases, 
variable biases and random errors. 

As presented in the body of the table, the trial consists of three 
engagements, each less than one minute in duration, four short bursts 
fired in rapid succession, and three rounds per burst. For Engagement 1, 
gun barrel surface temperature is at ambient air temperature and 
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therefore, barrel bend equals 0.0. The factor affecting azimuth fixed 
bias in this situation is the hard left to right cross-wind resulting in 
3.07 mils bias at the muzzle. Factors which affect the elevation fixed 
bias are the temperature differential from standard of 60°F to 0°F and 
the air density change from standard due to the cold temperature . The 
random error totals at the muzzle are the sum of the Round-to-Round , 
Burst- to-Burst , and Occasion-to-Occasion errors after use of the standard 
normal random number generator routine. Notice that the Occasion-to- 
Occasion errors do not vary within the engagement; the Burst- to-Burst 
errors vary only between bursts within the engagement; and the Round-to- 
Round errors vary with each round fired. Finally, the total miss dis- 
tance in meters at the target is expressed for each round in azimuth and 
elevation by considering both the system fixed bias and the random 
errors . 

In Engagement 2, the factors affecting the fixed biases of the wea- 
pon system remain the same; however, now barrel bend is introduced as an 
additional fixed bias . The 12 rounds fired in the first engagement 
affect the gun barrel surface temperature and the time between the two 
engagements, defined as the time for maximum bend to occur, has allowed 
the cool hard wind to decrease the outer surface temperature with the 
resultant bend. The new fixed bias in azimuth becomes 4.6?^ and 7 ■ 7^ 
for barrel bend and total fixed bias respectively. 

In Engagement 3» barrel bend is the only fixed bias change. The 
number of rounds fired in the first engagement (12) and the rounds 
fired in the second engagement (12) are cumulatively addressed to deter- 
mine the resulting outer surface temperature of the barrel. Although the 
same crosswind velocity is encountered, the effect of the higher outer 
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surface barrel temperature is a more severe barrel bend condition of 
7.01^. This results in a total fixed bias in azimuth of 10.08^1. 

There were a total of 480 trial runs conducted in the simulation 
model (24-0 trials with M60 fire control system approximation and 24-0 
trials with XM-1 fire control system approximation) . The results of the 
selected trial runs included in Appendix B adequately illustrate the 
effects of barrel bend in automatic cannon firing and validate the 
models developed. If additional trial run output is required for 
analysis, they are available upon request from the authors. 
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V. SUMMARY AND CONCLUSIONS 



Conventional warfare is a realistic possibility in the turbulent 
world of today . In a conventional war the tank will be a principle 
threat to allied forces and the primary counter to that threat. If the 
tank is to be effective in the antitank role, in view of the numerical 
superiority of opposing forces, there must be a capability of accurate 
fire well above that demonstrated in the Arab-Israeli War. Every attempt 
must be made to eliminate or minimize factors that affect the accuracy 
of tank gun fire . 

This paper has presented an analysis of factors affecting miss 
distance and proposes that barrel bend has been overlooked as a major 
contributor to inaccurate fire . A formula was derived for computing 
the expected value of bend which can realistically be experienced in 
varying weather conditions. A model was presented which computes barrel 
bend and other errors giving the resulting miss distance at the target. 

Analysis of the model output as presented in Appendix 3 provides 
the following conclusions. 

1. Miss distance at the target is less with the XM-1 fire control system 
approximation than the M60 fire control system approximation due primar- 
ily to the capability of sensors for air temperature, wind, and air 
density as well as the laser rangefinder. This minimizes significantly 
the system fixed biases other than barrel bend. The M60 fire control 
system is particularly susceptible to error induced by crosswinds. 

The fixed bias unaccounted for by the XM-1 fire control was generally 
less than 0.1 mil. This error is relatively insignificant when compared 
to the error from barrel bend. 
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2. The APFSDS round is more accurate than the HE round as a result of 
the projectile's design characteristics, e.g. significantly higher muzzle 
velocity . 

3. Rapid fire cannons -will be particularly susceptible to barrel bend 
due to the high barrel surface temperatures resulting from burst fire . 

4. Shorter bursts generate less heat and less barrel bend. 

5 . Cooler air temperatures increase the magnitude of barrel bend. 

6. Shorter barrels and barrels with thicker walls will be less sus- 
ceptible to bend . 

7. Barrel bend occurs during the cooling period after an engagement. 

The magnitude of bend that can occur depends on the number of rounds 
fired during the engagement, weather, and gun characteristics. An en- 
gagement of 20 rounds, i.e. four bursts of five rounds fired in rapid 
succession, is shown to generate sufficient barrel temperature to cause 
as much as a 12 mil error under wind cooling or nearly 16 mils under 
rain cooling. This bend becomes a fixed bias affecting the accuracy 

of rounds fired in the next engagement. 

8. In a series of engagements there will be an additive effect of 
barrel bend if sufficient time is not allowed for the bend to diminish. 
Thus the third engagement can be affected by the bend induced by the 
second engagement plus a portion of bend induced by the first engagement. 

There are ways to compensate for barrel bend and ways to eliminate 
bend altogether. If compensation for barrel bend is to be achieved, 
then extensive research is required to: 

1. Develop functions for outside barrel temperature. 

2. Develop functions for the cooling factor of varying amounts of wind 
and rain . 
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3* Determine whether the barrel does or does not straighten by any 
significant amount during firing, and if straightening does occur, at 
what cost in muzzle velocity from standard. 

4. Develop systems for sensing or direct measurement of barrel bend. 

On the other hand, if elimination of bend is desired, then research is 
required to develop a system which will insure even cooling. Shielding 
of the barrel from the elements of wind and rain by a protective shroud 
or coating seems to be the most viable alternative . This method is re- 
portedly being investigated by West Germany and Great Britain. 

In addition to further analysis of barrel bend, the testing of the 
models in the parametric analysis presented further areas for extended 
study . 

1. Examination of the effects of trunnion forces on muzzle position to 
provide additional input into the miss distance model. 

2. Determine the change in muzzle velocity from standard as a result 

of changes in propellant temperature due to heat transfer from the 

barrel temperature . 

3* Determine the change in muzzle velocity from standard as a result 

of the barrel bend condition and projectile straightening upon firing. 

This paper has particularly addressed the bend that can occur in 
rapid fire cannons because of the higher barrel surface temperatures 
that are likely to be experienced as these weapons are developed. The 
paper also concentrates on tank guns due to the stringent requirement 
for accuracy when firing at point targets with direct fire. However, 
bend occurs in weapons with lower rates of fire, including current tank 
guns and conventional artillery weapons. 

It is concluded that barrel bend can significantly decrease the 
accuracy of gunfire . The magnitude of the resultant error at the target 
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can no longer be tolerated if the tank is to remain the principal 
anti-tank weapon on the modem battlefield. Weapons under development, 
particularly rapid fire tank guns, should be required to demonstrate 
the capability of avoiding or at least compensating for barrel bend. 
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APPENDIX A 



ASSUMED FIRING TABLE AND EFFECTS DATA 

The firing table data presented in Appendix A is not the exact data 
being used for the ARPA 7 5 1 ™ cannon and selected ammunition, but is 
merely an assumed approximation for use in the models generated in this 
paper. For more effective evaluation of this developmental system, the 
exact firing table data should be used. 
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ASSUMED FIRING TABLE AND EFFECTS DATA 
ARPA APFSDS - 75MM 



RANGE SUPERELEVATION DRIFT 

(Meters) (MILS) (MILS) 

1500 2.8370 0 



2500 



4.8929 



0 



NON STANDARD - STANDARD 



RANGE 


CHANGE IN 


CHANGE IN 


NON-STD 


(M) 


HEIGHT (M) 


DRIFT (M) 


CONDITION 


1500 


.0060 


.0000 


1 m/s Velocity 




-.0025 


.0000 


1 Pet Density 




-.0015 


.0000 


1 Pet Air Temp 




.0000 


.0000 


1 m/s Range wind 




.0000 


.0358 


1 m/s Crosswind 




1 .4730 


.0000 


1 MIL Elevation 


2500 


.0144 


.0000 


1 m/s Velocity 




-.0110 


.0000 


1 Pet Density 




-.0030 


.0000 


1 Pet Air Temp 




.0015 


.0000 


1 m/s Range wind 




.0000 


.0983 


1 m/s Crosswind 




2.4540 


.0000 


1 MIL Elevation 
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ASSUMED FIRING TABLE AND EFFECTS DATA 
ARPA HE AMMO - 75MM 



RANGE 

(Meters) 

1500 

2500 



SUPERELEVATION 

(MILS) 

6 .6336 
12.8519 

NON STANDARD - STANDARD 



RANGE 


CHANGE IN 


CHANGE IN 


(M) 


HEIGHT (M) 


DRIFT (M) 


1500 


.0198 


.0000 




-.0218 


.0000 




-.0080 


.0000 




.0033 


.0000 




.0000 


*1970 




1 .4738 


.0000 


2500 


.0633 


.0000 




-.1298 


.0000 




-.0353 


.0000 




.0153 


.0000 




.0000 


.6140 




2.4553 


.0000 



DRIFT 

(MILS) 

0 

0 



NON-STD 

CONDITION 

1 m/s Velocity 
1 Pet Density 
1 Pet Air Temp 
1 m/s Rangewind 
1 m/s Crosswind 
1 MIL Elevation 

1 m/s Velocity 
1 Pet Density 
1 Pet Air Temp 
1 m/s Rangewind 
1 m/s Crosswind 
1 MIL Elevation 
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APPENDIX B 



OUTPUT TABLES 

The computer trial output tables included in this appendix represent 
a selected sample of the model output for analysis. The effect of each 
of the four selected weather conditions, which are constant for each of 
the eight tables, are examined by varying the trial conditions of the 
round, range, fire control and rounds in the burst. 
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COMPUTER TRIAL OUTPUT TABLE 
FIXED BIAS (MILS) 
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COMPUTER TRIAL OUTPUT TABLE 
FIXED BIAS (MILS) 
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COMPUTER TRIAL OUTPUT TABLE 
FIXED BIAS (MILS) 
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COMPUTER TRIAL OUTPUT TABLE 
FIXED BIAS (MILS) 
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APPENDIX C 



COMPUTER PROGRAM 

The computer program included in this appendix is that written for 
simulation of the miss distance and barrel bend models. The only changes 
required to simulate the various fire control system approximations 
(M60 or XM-i) are to change the data inputs for percentage of residual 
sensing errors and for Round- to-Round, Burst- to-Burst, and Occasion- to- 
Occastion dispersion. The program may also be adapted to different gun 
barrels by inserting the physical characteristics of the gun and an 
approximation for gun barrel surface temperature . 
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3R T PX I A I . R TR Y ( 4 ) »RTBX I A ) * QT0Y1 4) , OT 0 X ( 4 ) ,OTOY 14 ) , 
4CLMHTU).NRC(3) ,NPDSI 3) t A(6) f 3(23) ,C(20) 

DIFENSICN IFTU ,10) , t\< c T ( 15, 10) 

DC 10 J= l , l 5 

CL^ZILT-S^r OO . ) ).,A- X J/.P-4^^CWND(^-) r-RMNI.J U V AP RF-S U T^F P-i4 ) rFF (-2^3— 

1000 FCFMATC6F10.0) 

10 CCNT INLE 

or zc * = i 

P c A P ( E » 1 1 00 ) FTTPRO(K) r F TR NOE I K ) ,FT VLTY ( K ) , F TPMY (K ) f FTATMP (K ) , 
IFTRC-WC K) t F T C R W RI K ) , FT^LEVIK) , F T[)P F T (K) , F T SP E L I K ) , STOV c L(K) , 

2 A L 1 0 E F ( K ) 

1 100 FCFMA T ( I2F6 .0 ) 
cn M L T NU4- 



OC 30 1=1,4 

R C AC (5, 1200 ) RTPX(L),RTPY(L),RTBX(L) ,BTRY( L ) » C T f X ( L ) ,CTCy (L ) 
1200 PDF* AT ( 6 F 1 0 • 0 ) 

30 CCNTINLE 

P E *C ( 5 , 1300 )_ALEN, ALPHA,RD,R I , ( CUMH T ( I ) ,1=1,2) 

130 3 FCF M AT ( 6 F l 0 • 0 ) 

PFADI E, 1400 ) S TONTY, S T A Tm p r \C r RV<P t XCNAIM t XTNFCS»YTNFCS,V^L p T , 



85 
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Y3 = RT D Y ( I ) *C ( 2*M 
S l CX = X 1 +X2 + X3 

S I CY = Y 1 + Y2 * Y3 * 
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